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Abstract When a binary fluid demixes under a slow 
temperature ramp, nucleation, coarsening and sedimen- 
tation of droplets lead to an oscillatory evolution of 
the phase separating system. The advection of the sed- 
imenting droplets is found to be chaotic. The flow is 
driven by density differences between the two phases. 
Here, we show how image processing can be combined 
with particle tracking to resolve droplet size and veloc- 
ity simultaneously. Droplets are used as tracer particles, 
and the sedimentation velocity is determined. Taking 
these effects into account, droplets with radii in the 
range of 4 - 40/im are detected and tracked. Based on 
this data we resolve the oscillations in the droplet size 
distribution which are coupled to the convective flow. 



Keywords particle tracking • image processing • 
^ droplet size distribution • phase separation • binary 
^j" mixture 
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1 Introduction 

The characterization of particle distributions in fluids 
is important to control and optimize processes in food, 
pharmaceutical, oil and chemical industry, as Heffels 



|et al (19981) point out. To determine the particle mass 
flux, particle sizes and velocities have to be measured 
simultaneously, which was achieved by Petrak (2002). 



Here, we present a particle tracking algorithm, which 
uses droplets as marker particles to measure the flow 
field. They are created naturally in the phase separa- 
tion process of demixing binary systems. Droplet posi- 
tions and radii are detected simultaneously. The radius 
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can therefore be used as a criterion to identify droplets 
in subsequent images. Assuming Stokes law the sedi- 
mentation velocity can be calculated from the droplet 
radius. The droplet velocity is decomposed into sedi- 
mentation velocity and advection by the flow. By sub- 
tracting the sedimentation velocity from the droplet 
velocity, the advection of all droplets can be used to 
measure the flow field. With the advective flow field 
and the sedimentation velocity of each droplet, its po- 
sition in the next frame can be predicted and compared 
to the image. 

Previous studies have extensively used particle track- 
ing velocimetry as a measurement technique to inves- 



tigate turbulent flows carrying small particles (Maas 



et al 



(1993); 



Kreizer et al 



Malik et al (1993); Ouellette et al (2006); 



(2010)). To that end monodisperse tracer 



particles are added to the fluid whose position is de- 
tected by image processing of high speed camera data. 
While this procedure is followed successfully for mono- 
disperse particles in single-phase flows, it is not easily 
applicable for measuring flow patterns in demixing bi- 
nary systems, which have been studied by [Vollmer et al 



(1997 2002) and 



Emmanuel and Berkowitz 



(2006) 



Tracer particles cannot be added as they would act 
as nucleation centers, and therefore affect the droplet 
number density. Further more, colloidal particles aggre- 
gate on the interfaces and change the growth and coa- 
lescence rate, as |Thijssen and Cleg g (2010) have shown. 
For a review of the stabilizing effect of colloidal parti- 



(2003). 



cles in emulsions see Binks (2002) and Aveyard et al 



There also is a broad range of acoustic / electro- 
acoustic and optical techniques for the measurement of 



droplet size distributions, an overview given by |Maass 



et al (2009). Several laser based techniques measure 



chord lengths that have to be transfered into droplet 
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size distributions ( Hu et al] ([2006 ) ) , e.g. focus beam re- 
flectance measurement (FBRM, see e.g.|Ruf et al (2000)) 



or optical reflectance measurement (ORM, see e.g. 



et al 



et al 



2005| )). However, MaaB et al (2011) and Simmons 
(J2000|) have shown that these techniques give inac 



Lovick 



curate results for liquid/liquid dispersions in compari- 
son to image analysis of in situ microscope imaging. 

To investigate the evolution of demixing systems, 
it is therefore preferential to use particle tracking of 
droplets present in the demixing system, as described 
in the following. 

The article is organized as follows: We first char- 
acterize the chemical system. Then we introduce the 
experimental setup with a focus on the illumination 
technique, and explain the experimental procedure (sec. 
2). In section 3, we present the particle tracking algo- 
rithm. It is based on an image processing method to 
detect the droplets. In section 4, we measure the sedi- 
mentation velocity of droplets and give an indication on 
the uncertainties in the radius detection. We then use 
the velocity information to filter out falsely detected 
droplets. We show how the droplet size distribution and 
the characteristics of the flow field evolve in time. Fi- 
nally, in section 5, we discuss the applicability and the 
limitations of the droplet tracking procedure for the in- 
vestigation of demixing binary fluids. 



2 Experimental Method 

2.1 Model System 

As a model system a mixture of water and isobutoxy- 
ethanol (i-BE) is chosen. It has a lower critical point at 
25.5° C and demixes under heating. A phase diagram 
is shown in Fig. [I] The sample can easily be prepared 
at room temperature and demixes at temperatures be- 
tween 25.5°C and 80°C, as shown by |Nakayama et al 
( [2001] ). The i-BE (purity > 97%) is purchased from 
Wako Chemicals GmbH and used without purification. 
The data points in the phase diagram (Fig. HI) were 
determined by turbidity measurements. The solid lines 
show fits with sixth order (left branch) and fourth order 
(right branch) polynomials respectively. 

The density and the viscosity of the mixture depend 
on composition as well as on temperature. For the den- 
sity, the data of |Doi et al| ( |2000[ ) is used, taking thermal 
expansion and molar excess volume explicitly into ac- 
count. The determination of the viscosity is described 
in the appendix. 

Droplet detection is improved by using Nile Red as 
a fluorescent dye. Nile Red solved in butoxyethanol ab- 
sorbs in the wavelength range of 500nm to 590nm (green 
light) and emits between 580nm and 700nm (red). The 





50 


o 

o 


45 


CD 

~5 


40 


■i— > 




CO 




0) 


35 


Q_ 




F 




Q) 


30 



25 



\ w / 

\ i *\ /a 

X 2o>„ \y 


v -^„~-^ \ 



0.1 



0.2 0.3 0.4 0.5 0.6 
mass fraction $ of i-BE 



Fig. 1 Phase diagram of water/isobutoxyethanol mixture 
between 25 and 50° C. The samples are filled with the compo- 
sition ^ g . The mixture separates above a critical temperature 
T c = 25.5° C. For slow temperature ramps the compositions 
^b , ^t of the bottom and top phase follow the left and right 
branch of the binodal, respectively. 2<Po = <P t — ^b denotes the 
width of the miscibility gap. The left and the right branch are 
fitted with sixth and fourth order polynomials, respectively. 



dye solves preferentially in organic compounds and poor- 



ly in water. Fowler and Greenspan ( 1985 ) point out that 



the fluorescence of Nile Red depends strongly on the 
polarity of the solvent. In particlar, its fluorescence is 
quenched in water . Therefore, the i-BE-rich phase ap- 
pears bright and the water-rich phase stays relatively 
dark. Since each phase consists of a mixture of water 
and i-BE, the contrast between dark and bright parts 
in the images recorded during a measurement depends 
on the composition of the two phases. It was checked 
that the dye has no significant influence on the phase 
diagram. 



2.2 Experimental Setup 

The probe is contained in a fluorescence cell (10 x 10 x 
33mm) 117.100F-QS made by Hellma GmbH. The mea- 
surement cell is mounted in a water bath with con- 
trolled temperature (see Fig. [2|. An immersion cooler 
Haake EK20 is cooling with constant power, and a tem- 
perature control module Haake DC30 is heating the 
water bath to a preset temperature. Additionally the 
temperature of the water near the sample is measured 
with a PT100 temperature sensor. The temperature is 
controlled with an accuracy of 15 mK. 

To get enough signal from the fluorescently labeled 
i-BE-rich phase, a bright light source is required. A 
laser light sheet is not suitable for illumination since 
it is unidirectional and the two phases are not index- 
matched. Each droplet acts as a little lens which fo- 
cuses and diffracts the parallel laser light sheet. After 
passing a short distance of the sample, the droplets 
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have transformed the uniform light sheet into a stripe 
pattern, which illuminates the droplets very inhomoge- 
neously. In contrast a mercury short arc lamp (LOT- 
Oriel, 100W) with a bright light emitting spot gives far 
better results. It provides a high light intensity, which 
is diffuse enough that almost no stripes appear (figure 

IS- 

A sketch of the optical components is given in fig- 
ure [2] The light is collected by the collimator lens (C) 
and a bandpass filter (GF) selects the green emission 
lines (546 and 577/579nm) for the excitation of the 
fluorescent dye. The spherical lens (LI, / = 200mm) 
forms a parallel light beam. The cylindrical lens (L2, 
/ = 80mm) focuses the light normal to the focal plane 
of the camera. A spherical diaphragm (D) allows to ad- 
just the amount of light manually. With a narrow slit 
of 300/im width directly in front of the measurement 
cell a vertical light sheet is formed. The measurement 
cell is covered by black apertures to shield stray light. 
The fluorescent light from the illuminated plane of the 
sample is projected by a / = 35mm objective (L3) to 
the chip of a BM-500CL monochrome progressive scan 
CCD camera. It takes 2058 x 2448 pixel images with a 
maximum frame rate of 15Hz, covering 1. 3mm x 1.5mm 
of the sample (1.6 pixel//xm). With a longpass filter 
(RF, edge 594nm) the excitation light is filtered. 

The camera, temperature control and magnetic stir- 
rer are connected to a computer. Measurements are 
fully automated using a LAB VIEW program. 



2.3 Experimental Procedure 

The temperature protocol for our experimental proce- 



dure is adapted from Auernhammer et al ( 2005 ) . After 



mixing the sample at 24° C for one hour, the tempera- 
ture is set to T = 25.8° C (0.3K above critical point) and 
the system equilibrates for four hours. Then a tempera- 
ture ramp is run from 25.8°C to 50° C, which is designed 
to keep the volume flux between the phases constant. 



Following Auernhammer et al (2005) and Cates et al 



( 2003 ) a driving parameter £ is defined 



f = 



1 



dt 



(1) 



*o(T(t)) 

with 2^o = ^t — ^b? where $ t an d ^>b denote the com- 
positions of the top and the bottom phase, respectively 
(compare Fig. IT]). Hence, 2<£>o(T) is the width of the 
miscibility gap for a given temperature T. The driving 
parameter £ is equal to the rate of droplet mass pro- 
duction in the two phases. Equation ([!]) is inverted to 
calculate a temperature ramp with fixed driving, 

£*o(T(t)) 



measurement cell 

water 
bath 




T(t + Si) = T(t) 



&r$o(T(t)) 



St. 



(2) 



camera 



Fig. 2 Sketch of the experimental setup, top view. The light 
of a short arc Hg vapor lamp is collected by a collimator lens 

(C) and parallelized by a spherical lens (LI). A cylindrical 
lens (L2) forms a light sheet which is cut by a diaphragm 

(D) and a line aperture (300//m width) directly in front of 
the sample cell. The light sheet illuminates a plane of the 
binary mixture, a green filter (GF) selects the wavelength 
range suitable for excitation of the dye. The emission of the 
fluorescently labeled phase is passes a red filter (RF) and is 
projected by an objective (L3) onto a 5 Mpixel CCD cam- 
era chip. The sample cell is mounted in a water bath with 
controlled temperature. 



Equation ^ is integrated with 10000 uniformly dis- 
tributed time steps in the range [to,t en d] chosen such 
that T(t ) = 25.8°C and T(t end ) = 50.0°C. To im- 
plement the temperature scan, the temperature of the 
thermostat is set with a rate of 1 Hz according to the 
calculated temperature ramp. 



3 Particle Detection and Tracking 

In this section we describe a particle detection and 
tracking algorithm for polydisperse sedimenting parti- 
cles. Radii and positions are detected simultaneously. 
The first part describes an image processing algorithm 
which detects droplets in separate images. In the next 
step the flow field is estimated. Then the droplets are 
tracked through a series of images. Finally the velocity 
field is recalculated based on the droplet trajectories. 
This procedure allows us to sort out artifacts (caused 
e.g. by dirt or overlaps of droplet images) of the image 
processing, and to determine Lagrangian particle veloc- 
ities. The image processing is done with the MATLAB 
Image Processing Toolbox. The data (typically 20000 
images per measurement) is processed on a computer 
cluster. For simplicity we restrict our description to the 
detection of the fluorescently labeled droplets in the 
bottom phase. By inverting the images and the direc- 
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tion of gravity the same algorithm can be used for the 
dark droplets detected in the top phase. 



3.1 Processing of Single Images 

3.1.1 Preprocessing 

In a first step, dark spots caused by dirt on the camera 
chip are removed with a flat field correction. Next the 
image is Fourier filtered and the contrast is optimized. 
To reduce the horizontal stripes (illumination from the 
right) produced by the light dispersion of the droplets 
a high pass filter (Gaussian filter with width 5 in hori- 
zontal and 200 in vertical direction) is applied. For re- 
ducing the noise, a low pass filter (isotropic Gaussian 
filter, a = 100) is applied. In figure [31 a raw image and 
a preprocessed image of fluorescent i-BE-rich droplets 
in the denser water-rich bottom phase are shown. For 
high droplet densities, there is a significant number of 
overlaps of droplet images. Therefore we need special 
algorithms to take care of this for detection. 

3.1.2 Droplet Detection 

The droplets are detected in the preprocessed image 
using two strategies. The first relies on thresholding 
the image, the other uses marker-controlled watershed 
segmentation. The first strategy is visualized in figure 
[4] giving an example. All pixels darker than a thresh- 
old level (Fig. 4(a) ) are identified as background (Fig. 
|4(b)| ). Isolated background regions are deleted. Then, 
for each pixel the distance to the background is cal- 
culated (Fig. |4(c)[ ). Local maxima in this distance map 
correspond to droplet centers, their distance to the back- 
ground being their radius (Fig. |4(d)| ). This procedure is 
carried out twice with two thresholds, one being better 
for finding big droplets, the other for small droplets. 

The second strategy is based on the marker-control- 
led watershed segmentation algorithm proposed by|Gon-| 



zalez et al 



(2004) pp. 422. It is illustrated in figure pM 



Extended maxima in the image (Fig. |5(a)| ) are used as 
marker for the droplets (Fig. |5(b)[ ). Lines between these 
maxima as well as extended minima are markers for the 
background. With a sobel filter the intensity gradient 
in the image is calculated (Fig. |5(c)"| ). The markers are 
imposed as minima on the gradient image upon which a 
watershed transform is operated. The watershed lines of 
the intensity gradient correspond to droplet edges. The 
markers are needed to avoid over-segmentation. Areas 
and centroids of the regions enclosed by the watershed 
lines are taken for the droplet radii and positions (Fig. 




Fig. 3 (a) Raw image and (b) preprocessed image of i-BE- 
rich droplets in the water-rich bottom phase with f = 1.05 • 
10~ 5 s at 25.98 ± 0.02°C. Horizontal stripes and the vertical 
intensity gradient are removed by Fourier filtering. 



Using these two approaches, even small droplets on 
top of bigger ones and overlapping images of droplets 
can be resolved. On the other hand, some droplets are 
found multiple times and some imaging artifacts are 
taken for droplets. If two droplets with radii r\^r^ are 
detected close to each other (centroid distance less than 
one radius) with almost the same radius (ratio 0.7 < 
7"iA*2 < 1-4), they are considered as two detections of 
the same physical droplet. 
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Fig. 4 Sketch of distance based image segmentation. 

(a) False color plot of two overlapping images of droplets. 

(b) Thresholded image, (c) Distance from background. 

(d) Detected droplets (black) based on local maxima of dis- 
tance. 



3.1.3 Matching to Image 

Not fully compensated dirt on the camera chip, overlaps 
of droplet images or unsharp droplet images can cause 
false detections. Hence, a second step is implemented to 
distinguish droplets from artifacts and the estimates of 
the droplet radius r and position (x, y) [pixels] are com- 
pared to the preprocessed image. As characteristics of 
droplets in the image, we consider them being brighter 
than their environment and having a circular intensity 
gradient at the edge (figure [6|. A match function m is 
defined as 

m(x, y, r) = A(r) [IdropO, V, r) - /ring (a, 2/, r)] 

B (x-x r ) 2 + (y-y r ) 2 ( 3 ) 

where the radius r and the coordinates x, x ri y and y r 
are given in pixels. The difference between the mean 
intensity of the droplet /drop {x, y,r) and the mean in- 
tensity J r i ng (x,2/,r) of a ring of 3 pixels width around 
it measures the intensity gradient at the droplet edge 



(figure 6(a)). An empirical radius dependent prefac- 
tor A(r) = 1/ (30 + r 13 ) ensures comparable matches 
for small and big droplets. Often small droplets are 
detected erroneously at the edge of big droplets. In 
this case, a little part of its halo is covering the dark 
background and the remaining part is covering the big 
droplet. The intensity weighted centroid of the halo (de- 
noted as x r and y r ) is calculated and compared to the 



Fig. 5 Sketch of watershed image segmentation, (a) False 
color plot of two overlapping images of droplets, (b) Bright 
extended maxima mark droplets, lines in between and dark 
regions mark background, (c) Markers are imposed as minima 
(blue) on gradient image. Watershed lines follow droplet edges 
(red) . (d) Detected droplets (black) based on regions enclosed 
by watershed lines. 





Fig. 6 Definition sketch for the match function, (a) /drop 
denotes the mean intensity of the droplet, / r in g the intensity 
of a ring around it. (b) (x,y) is the centroid of the droplet 
with radius r, (x T ,y r ) is the intensity weighted centroid of the 
ring. 



droplet position (figure |6(b)[ ) . The deviation accounts 
for this asymmetry and reduces the value of the match. 
It is weighted with the empirical prefactor B = 3. For 
every droplet found in the droplet detection step, the 
match with the image is calculated. By varying radius 
and position of the droplet a local maximum of the 
match is found. All droplets with a match bigger than 
0.08 are used for further analysis. 

A final step to increase the reliability of the mea- 
sured size distribution takes advantage of the time in- 
formation of image sequences. Tracking the droplets al- 
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lows to sort out unphysical trajectories of artifacts and 
to find undetected droplets in the image using its past 
or future trajectory. To implement this information in 
section 13.3.31 we must evaluate the flow field and the 
droplet trajectories in the flow. 



3.2 Calculation of Flow Field 

3.2.1 Identification of Droplets in Consecutive Images 

Typically droplets with radii 4/im < r < 40/im are de- 
tected. While small droplets (r < 10/im) closely follow 
the flow, large droplets (r > 20/im) are mainly driven 
by their buoyancy and do not qualify as tracers. They 
sediment towards the interface due to the density dif- 
ference of the two phases. In a Lagrangian frame, co- 
moving with the surrounding fluid, the sedimentation 
velocity u se( ± amounts to the Stokes velocity for a sphere 
of radius r, which is (Batchelor, 2005, p. 234) 



^sed 



2 Ap gr 2 
9 rj ' 



(4) 



Here Ap denotes the mass density difference between 
sphere and fluid, g the gravitational acceleration, and 
r] the dynamic viscosity of the fluid. Hence, the droplet 
velocity i^drop can be split into a sedimentation and an 
advection term 



^dropO, V, r, i) = u sed (r) + iiflowO, y, t). 



(5) 



This ansatz will be checked in section 14.11 To track 
droplets the corresponding images of droplets in con- 
secutive frames must be identified. To this end, the ve- 
locity field of the previous time step Uf[ ow (x,y,t — St) 
is taken as an initial guess for u^ ow (x 1 y 1 t). For each 
droplet found in one image the position in the next im- 
age is predicted using equation ([5|. The prediction is 
compared to the droplets found in the image, and the 
droplet pair with the closest distance between predic- 
tion and actual position and with similar radii is iden- 
tified as one physical droplet. In figure [7] the droplets 
in the image are marked with red circles and the corre- 
sponding droplets in the next image with green circles. 
The procedure used to determine the flow field, which 



is indicated by the arrows, is lined out in section 3.2.2 



Subsequently in section |3.3| we describe details of the 
scheme adopted for finding particle trajectories. 

3.2.2 Determine the Flow Field 

In order to calculate the Eulerian flow field Ufi ow (x, y, t) 
the image is sampled by an equidistant mesh with 25 x 
20 grid points. The grid distance is about 60/im, which 
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Fig. 7 Fl uorescently labeled droplets in the bottom phase 
(Fig. |3(bj] ). The arrows correspond to the flow field. Red cir- 
cles mark the droplets found in the image and green circles 
their position and radius in the next image. 



is in the order of the droplet distance. The correlation 
length of the flow field is about 200/im, so the grid 
can resolve the flow structures. A flow correlation time 
of 30s allows for small frame rates (order of 1Hz). By 
subtracting the sedimentation velocity u se( ± from the 
droplet velocity i^drop the flow contribution Wfl ow can 
be calculated: The advection of the N droplets in each 
mesh cell (i) is averaged, incorporating with half weight 
the M droplets of neighboring cells (ij): 



x flow 



1 



7V + M/2 



N 

E( 



drop 



(0 

i sed 



fc=l 



M 4 

k=ij=i 



(6) 



k (ij) _k AiJh 
"-drop fl 



^sed, 



3.2.3 Smoothing of Flow Field 



.« 



The advection field u^ w is smoothed with a weighted 
average of next-nearest neighboring cells to eliminate 
unphysical discontinuities. 



x flow 



1 



x flow 



1 



STvpl 



1 



i=i 



flow 



— V u (ij,) 



%w 



(7) 



next neighbours next-nearest neighbours 



In figure [7| a smoothed advection field is shown together 
with the matched droplets. The Eulerian flow field is 
used to predict the droplet positions in the next image. 
Note that the spatial distribution of detected droplets 
is not homogeneous and some grid points can only be 
calculated by interpolation. For the subsequent analysis 



Particle tracking for polydisperse droplets 



of the flow field only the grid points which are directly 
computed from droplet displacements are taken into ac- 
count. 



3.3 Particle Tracking 

3.3.1 Forward Tracking 

Now that the advection field and sedimentation veloci- 
ties for each droplet found in the image at time t have 
been determined, the position at time t + St can be 
predicted. This prediction is matched to the image of 
t + St and a local maximum of the match is searched 
for. This procedure is illustrated in figure [8j Having de- 
termined the displacement from time t to t + St, the 
advection field can be recalculated. In the next step, 
the positions of droplets tracked from image t to t + St 
plus all other droplets also detected in image t-\-St are 
predicted for the time step t + 2St. They are compared 
to the droplets in image t + 2St and the procedure is 
repeated. The identified droplet pairs are sorted into 
trajectories, which contain the position and radius val- 
ues for the different frames. Based on the trajectories, 
the flow field Uf[ ow (x,y,t) is recalculated according to 



the procedure described in section 3.2.2 





Fig. 8 Forward tracking: (a) Sedimentation velocity (blue, 
upward) and advection (green, 2 o'clock), (b) predict the next 
position (red dashed), (c) match the prediction to the next 
image (red solid) and (d) build trajectory (black). 



3.3.2 Backward Tracking 

The forward tracking result can be improved by subse- 
quently tracking the droplets backwards in time. The 
positions of droplets at the beginning of the trajecto- 
ries are predicted for the time step before and matched 
to the image (compare figure [9]). The droplet is com- 
pared to the entries in the other trajectories to merge 
interrupted trajectories, and to suppress branching of 
trajectories. As yet, our tracking algorithm does not 
account for the collision of droplets. 





Fig. 9 Backward tracking: (a) Sedimentation velocity 
(blue, upward) and advection (green, 8 o'clock), (b) predict 
the previous position (red dashed), (c) match the prediction 
to the previous image (red solid) and (d) enlarge the trajec- 
tory (black). 



3.3.3 Filtering of Trajectories 

Trajectories shorter than three time steps are removed. 
Finally, all droplets which have moved by less than half 
the average advection velocity during the tracked time 
interval are deleted. Most likely, those spots are imaging 
artifacts. Eventually, the flow field is recalculated once 
more based on the enlarged and filtered trajectories in 
order to arrive at a consistent data set. 
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4 Results 



4.1 Sedimentation Velocities 



We first aim to confirm the validity of the ansatz Q and 
([5]) . In figi: 



mre 



10 the vertical velocity calculated from the 
trajectories of the tracked droplets is shown as a func- 
tion of the radius (red data points). The error bars dis- 
play the uncertainty of the radius detection (see below) 
and the change of the material parameters in the time 
interval of averaging. It is compared to the sedimenta- 
tion velocity (blue points) due to equation [4] modified 
by the mean vertical velocity of the flow, which is cal- 
culated from the flow field. The Stokes velocity of a 
sphere Q appears to describe the average settling rate 
very well. We attribute the deviations for small droplet 
radii to the fact that both the droplet size distribu- 
tion and the average flow velocity oscillate. Most of the 
small droplets can be found at the beginning of the os- 
cillations, when the average downward flow velocity is 
relatively small (see below). 
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Fig. 10 Comparison of the measured vertical velocity of the 
droplets (red) and the Stokes velocity modified by the mean 
flow (blue), in the bottom phase at 26° C. 



The correlation time of the flow is found to be Tfl ow ~ 
30s. As an upper limit for the droplet response time 
we find Td rop = 36/is at T = 50° C for droplets with 
a = 40/im. Therefore, typical Stokes numbers in the 
experiment are smaller than 10 -6 . Even for a flow with 
high turbulence intensity with a Kolmogorov time scale 
of the order of microseconds, the Stokes number would 
not exceed 10 -3 . Due to the small density differences of 
fluids, kinematic particle-turbulence interactions (for a 
review see Vaillancourt and Yau (2000)) can hence be 



ruled out for droplets, even in a turbulent flow. 



4.3 Trajectory Length 

For a quantitative analysis trajectories of 2000 images 
are sorted into radius bins according to their average ra- 
dius. In each bin the number of trajectories of a given 
length is color coded on a decade logarithmic scale (fig- 
11). The red line indicates the time t = H/u se d 



ure 

needed by a droplet settling with Stokes velocity to pass 
through the image of height H. It puts an upper bound 
on the trajectory length and leads to a decrease of the 
trajectory length for big radii. 

For small radii, the duration of individual droplet 
detection is observed to increase roughly linearly with 



the radius (figure 11 ). Hence we assume that the depth, 



in which a droplet can be detected, is proportional to 
its radius. Small droplets have to be exactly in the 
focal plane of the camera to be detectable, where as 
big droplets can still be detected when their center is 
slightly off the focal plane. Since the radius is detected 
by the maximum of the intensity gradient, it is rather 
robust to slight blurring, and once the droplets are too 
far from the focal plane, they are no longer detected be- 
cause the intensity gradient at their edge is below the 
threshold of the match function. 



4.2 Stokes Number 

Equation ([5| holds for droplets with a Stokes number 
St <C 1. The Stokes number is defined as the ratio of 
two time scales 



St = 



7"drop 
7~flow 



(8) 



In a flow which changes on the time scale Tfl ow the ad- 
vected droplets respond to accelerations on the time 
scale 



7"drop 



2 Apr 2 
9r] 



(9) 



4.4 Uncertainty in Radius Detection 

The fluctuations of the detected droplet radii in the 
trajectories are calculated to estimate the uncertainty 
of the attributed radius. For radii r < 19/im a linear 
increase of the standard deviation is observed, yielding 
a constant relative uncertainty of about 20%. Above 
19/im, the standard deviation is more or less constant 
with 3.8/im. 



4.5 Oscillating Size Distribution of Droplets 

By slowly heating the system, more than twenty os- 
cillations in the size distribution can be found (Fig. 
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Fig. 11 Radius dependence of the trajectory length: The 
number of trajectories with a given radius and length are 
color coded on a logarithmic scale. The black squares denote 
the average length for each radius bin and the solid red line 
is the sedimentation boundary. 
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Fig. 12 Radius-dependent standard deviation of the radius 
of tracked droplets.. 
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Fig. 13 Size distribution of droplets of the bottom phase 
with £ = 1.05 • 10 _5 s _1 : (a) before and (b) after consistency 
check based on droplet tracking. The droplet number den- 
sity per radius given in /xm -4 is color coded on a decadic 
logarithmic scale. 



4.6 Oscillating Flow Properties 



13). A complex interplay of nucleation ( Krishnamurthy 



2008 



2007 



and Goldburg (1980); Sagui and Grant|( 1999); 



)), coarsening (Aarts et al ( 2005 ) ; [Vollmer et al 



)) and sedimentation (Hayase et al (2008)) causes 



Vollmer 



To characterize the time evolution of the flow field (us- 
ing only non- interpolated grid points), the mean verti- 
cal and horizontal flow velocity as well as the root mean 
squared velocity components are calculated. For aver- 
aging in time, four bins for each oscillation are used (fig- 



oscillations of the supersaturation and turbidity, mod- ure 14). The mean horizontal velocity fluctuates around 



eled by Vollmer et al (2007) and Benczik and Vollmer 



(2010). They can also clearly be identified in false color 



plots of the droplet size distribution displayed in Fig. 
[l3j Comparing the two size distribution plots in Fig. [13 
reveals the effectiveness of droplet tracking. The noise 
in the processing of single images is reduced signifi- 
cantly and the evolution of the size distribution is un- 
raveled for small droplets also. A detailed analysis of 
the oscillating size distribution and its implications for 
the understanding of precipitation in other systems (for 
example rain formation, for a review see Shaw (2003)) 
will be given elsewhere. 



zero as expected for symmetry reasons, and the vertical 
flow velocity is typically negative. This can be explained 
by two effects. On the one hand, the upward sediment- 
ing droplet volume has to be balanced by a downward 
fluid motion to ensure a zero net volume flux. On the 
other hand, there is a large scale convection since the 
system is heated from outside. Close to the walls the 
compositions of the phases are slightly more separated. 
Hence the density of the bulk is higher close to the wall 
than in the central region of the sample. This drives a 
convection pattern with downward motion close to the 
wall. A coupling of the oscillating droplet size distribu- 
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tion with the flow field produces oscillations in the aver- 
age vertical flow component and the root mean squared 
velocities of both components. 




time[h] 

Fig. 14 Mean fluid velocity in horizontal (red circles) and 
vertical (blue circles) direction, and the rms velocity in hor- 
izontal (red triangles) and vertical (blue triangles) direction. 
Data of twelve oscillations are shown, where four consecutive 
data points correspond to one oscillation. The measurement 
is performed in the bottom phase with £ = 1.05 • 10 _5 s _1 . 



4.7 Oscillation Periods 

In each measurement, the driving rate £ is kept approx- 
imately constant. One may therefore expect the system 
to react with a constant time scale. However, we find, 
that during one experimental run, the oscillation pe- 
riod At decreases with increasing temperature (com- 



pare figures 13 and 14). This trend has been observed 
in turbidity measurements with other binary mixtures 



by Auernhammer et al (2005) and Vollmer et al (2007) 



before. Secondly, the oscillation period depends on the 
driving rate. For faster temperature ramps, the oscilla- 
tion period decreases. We have measured the oscillation 
periods of the droplet size distribution for two decades 
in £ for the top phase and almost three decades in £ for 



the bottom phase (see figure 15). The colors and sym- 
bols encode different heating rates £: black stars, £ < 
6xl0" 6 s _1 ; blue crosses, £ < 1.3xl0 _5 s _1 ; cobalt cir- 
cles, £ < 3 x 10 -5 s _1 ; green triangles, £ < 6 x 10 -5 s _1 ; 
red squares, £ < 3 x 10 _4 s _1 ; and magenta diamonds, 
i > 3 x 10" 4 s _1 . A trend of the period At with the 
driving £ is visible, but it is masked by the tempera- 
ture dependence of At. 
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Fig. 15 Oscillation periods At as a function of the driv- 
ing rate £ for the top (a) and the bottom (b) phase. In 
both phases, the oscillation period is smaller for higher driv- 
ing rates. In each experimental run, the driving is kept ap- 
proximately constant. Nevertheless, the oscillation period de- 
creases with increasing temperature. The use of the symbols 
is described in the text. 



5 Discussion &; Conclusions 

To investigate the physics of phase separation we study 
the water/isobutoxyethanol-system. It has very conve- 
nient physical properties, a critical point just above 
room temperature, an upper miscibility gap, and low 
vapor pressure. Therefore, sample preparation is easy 
and a broad range of parameters can be studied with a 
fully automatized setup. Using Nile Red as a dye and 
choosing a system with one polar and one organic phase 
we obtained images with a high fluorescence contrast. 

A short arc mercury vapor lamp is suitable for the 
illumination of the sample. To follow the evolution of 
the broad droplet size distribution for several hours, 
high resolution images (5 Mpixel) are taken at small 
frame rates (0.3 to 5 Hz). The resolution of the camera 
has to be large enough to detect droplet radii, whereas 
the frame rate must be adapted to the maximal droplet 
velocities. For measuring well resolved flow fields the 
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droplet number density has to be high enough. This 
sets a lower bound to the rate of temperature change, 
which controls the nucleation of droplets. On the other 
hand, for high droplet number densities, measurements 
are only possible close to the walls. As the droplets are 
not index matched, clear images can only be obtained at 
a wall distance up to once or twice the mean distance of 
droplets. Therefore the measurements were conducted 
0.3mm away from the wall, i.e. at a distance of about 
1.5 times the correlation length calculated from the flow 
fields. 

We presented a particle tracking algorithm based on 
image processing with MATLAB. It detects droplets in 
a range of 4 to 40 jam radius (bigger ones can be de- 
tected but are very scarce) in a sample field of 1.3 x 
1.5mm 2 . The motion of the droplets is decomposed into 
a sedimentation and an advection term. The sedimenta- 
tion velocity is given by the Stokes velocity of a sphere. 
The advection of droplets is sampled on an Eulerian 
grid. By taking the sedimentation velocity explicitly 
into account, big droplets can also be used to calculate 
the flow field, which reduces the noise. This is possible 
because the Stokes number is much smaller than one 
for all droplets in our experiments. It is not necessary 
to know in advance the radius at which sedimentation 
starts to become significant. This threshold depends on 
the specific flow conditions. 

There are two advantages by incorporating informa- 
tion on the droplet radius: First, the radius is a good 
criterion to identify droplets in subsequent images. Sec- 
ondly, the prediction of the droplet position in the next 
image is improved by knowledge of the radius depen- 
dent sedimentation velocity. 

The experimental setup and the tracking algorithm 
are optimized for long measurement times (several hours) 
and low Reynolds number flow. Our experimental pro- 
cedure can be applied to a broad variety of binary mix- 
tures, provided an appropriate fluorescent dye for la- 
beling (only) one of the two phases is found and the 
excitation and emission filters are adapted to it. It en- 
ables then detailed investigations of the evolution of 
the system: The droplet size distribution, Lagrangian 
particle velocities and Eulerian flow fields can be mea- 
sured simultaneously for a broad range of heating rates, 
temperatures and sample geometries. 

To investigate other laminar or turbulent flows a 
high speed camera has to be used. Our technique might 
then be a promising approach to investigate turbulent 
boundary layers with reactions or phase separation. Fur- 
thermore, one may consider placing an endoscope into 



the sample, as Maass et al (2009) have done for exam- 



ple, to measure the size distribution outside the bound- 
ary layer. 



A further advantage of the technique is that the si- 
multaneous measurement of particle position and par- 
ticle size allows us to determine the size evolution of 
individual droplets as they progress across the mea- 
surement area. This information can be used to ob- 
tain information on coalescence rates and collision ef- 
ficiency of sediment ing droplets: A preliminary study 
shows that droplets with radii about 40/im grow on av- 
erage by 0.5/xm/s while they travel across the measure- 
ment area. This growth amounts to a collision efficiency 



(see Pruppacher and Klett (1997) pp. 569 and Pinsky 
et al (1999)) of order unity. Further studies are under 



way to determine how the collision efficiency depends 
on the radius of the sedimenting droplet. This method- 
ology should be of interest to precipitation processes 
such as rain formation, where the growth of medium 
sized droplets is believed to be driven by coalescence. 
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Appendix 

The viscosity of i-BE is measured with an Ubbelohde 
viscosimeter type 537 10/1 made by Schott. The tem- 
perature dependence of the viscosity r] [kg/ms] is fitted 

by 



V(T) 



B • (20- T)-C-(T- 20) Z 

A • 10 T + D 



(10) 



with temperatures expressed in °C. The coefficients are 
given in tabl e [l] where the values for water are taken 
from iWeastl (11988-1989). 



A[kg/ms] 



B 



C 



D 



water 
i-BE 



1.002- 10^" 
3.36- 10- 3 



1.3272 
1.730 



0.001053 
0.001 



105 
108 



Tab le 1 F it coefficients for the viscosity of water (Weast 
( |1988-1989| )) and i-BE, defined by equation {To} . 



To interpolate the viscosities for a mixed phase of 
given mass fraction <& we use the composition-dependent 
viscosities at 25° C for a homogeneous mixture in the 



single phase regime given in Menzel et al (2003). The 
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Fig. 16 From top to bottom: Viscosity of the top phase, 
isobutoxyethanol, the bottom phase and water, using equa- 
tions ( |10|12| ). At 40° C the viscosity of the two mixed phases 
was measured with an Ubbelohde viscosimeter. 



data is fitted with a fifth order polynomial 



/?(<£, T = 25°C) 



-4O.660 5 + 1O3.440 4 - 100.32<Z> 3 
+39.350 2 + 0.170 + 0.91. 

(11) 



Assuming that the coefficients of interpolation are 
not changing substantially in the temperature range of 
our measurements, a rescaled viscosity fj(@) is defined. 
It only depends on the composition <P 

V (*,T) = fj{$) ■ 1! % be{T) + (1 - fH$)) ■ VH 2 o(T)- (12) 



The viscosities of the two phases are shown in figure 16 
as a function of temperature. To check the strong as- 
sumption entering this interpolation, we measured the 
viscosity of the two phases at T = 40° C. For both 



phases it followed the prediction of equation (12) to 



within 2%. This is sufficiently accurate for our means. 
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